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SUMMARY 



The result of the present investigations Is substan- 
tially as follows. The buckling of centrally loaded col- 
umns of open section 1b always accompanied by a twist if 
the cross section discloses neither axial nor point symme- 
try. Then the cross sections twist about an axis of ro- 
tation D, the location of which depends upon the shape 
of the median line of the section, the wall thickness and 
column length, and the limiting conditions. There are 
three such axes and consequently three different critical ' 
compressive stresses (twisting failure stresses). With 
point symmetry of the cross section (wherein the case of 
double-column symmetry is contalnod as a special case) the 
three orltical compressive stresses are givon in two Eulor 
strossos for (twist-free) buckling in dlroction of the two 
principal axos of inertia and one twisting failure stress 
for twisting about an axis of rotation passing through the 
contor of gravity. With simple cross-section symmetry, it 
finally affords one Euler stress for buckling in direction 
of the axis of symmetry and two twisting failure stresses 
for twisting about two axes of rotation in the plane of 
symmetry. Ruckling perpendicular to the axis of symmetry 
is therefore connected with a twist of the column. The 
thicker the wall and the greater the length of the columns 
the more the effect of the twist is neutralised, as the 
distance betwoen center of rotation and center of gravity 
continues to increase until finally, the observed buck- 
ling is practically free from twist. 

This holds for symmetrical as well as for unsymmetri- 
cal sections; the Euler formula gives. In this case, good 
(slightly too high) approximate values. 

•"Drillknicken zentrisch gedruckter Stabe mit offenen Pro- 
fil. im elastischen Sereich." Luf tf ahrtf or schung, vol. 
14, no, 9, September 20, 1937, pp. 444-457. ■ 
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INTRODUCTION 



For centrally loaded columns of solid cross section 
or thick-wall ed, hollow cross section, only one type of 
instability of the "straight" equilibrium pattern is rec- 
ognized r namely, twist-free buckling, which is usually 
termed flexural buckling or (in the elastic range) Suler 
buckling. For cross sections of thin-walled open or 
cloned soction and short column length, another type of 
instability in the original equilibrium pattern is of im- 
portance, namely, the buckling of the less bending-resis- 
tant walls due to inferior wall thickness. Lastly, there 
is yet a third instability phenomenon observed on thin- 
walled columns of open section - that is, torsional buc k- 
Una or, aB it is also called, twi sting failur e. Here, 




under a certain critical compressive force, the straight 
equilibrium pattern may be accompanied by infinitely ad- 
jacent twisted equilibrium patterns whereby, in contrast 
to tho^ case of buckling, tho cross-sectional shape is pro- 
served" and the strains due to reciprocal twisting of the 
cross sections about a well-dof inod axis of rotation are 
given. 

The roason for the importance of twisting failure be- 
ing rostrictod to columns of thin-walled open section, is 
duo to the fact that opon sections - especially, when 
thin— walled — possess an oxtremoly low twisting strength. 
Twisting failuro is most often and most clearly observed 
in airplane designs where the employed sections usually 
have nuch thinner walls than is customary on other struc- 
tures. Aviation literature on this stability problem is 
substantially exhausted with a theoretical treatiso by H. 
Wagner (reference l), and a further report by H. Wagner 
and T7. Pretschner (reference 2) which also contains test 
data for angle sections. In structural engineering lit- 
erature, the twisting failure of centrally loaded columns 
doos not appoar at all - as far as the writer known - ex- 
cept for one recent article by H. and F. Bleich (referonco 
3). In this report tho writors use tho energy method for 
the derivation of the differential equations for the strain 
quantities, while Wagner attains a clearer differential 
equation for the angle of twist from a consideration of 
the moment equilibrium about the column axis. In spite of 
the fundamentally identical assumptions, the results do 
not agree, and for the following reasons: 

Wagner defines the axis of rotation, to begin with, 
on the inconclusive assumption that the center of rotation 
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on twisting, fail-are is coincident with the so-called "cen- 
ter of -shear ," while--H. and 3!« Bloich 'a formula for the 
work of the. external forces, fails to include the term 
which accounts for the twisting of the cross sections. 

In the following a complete theory of twisting fail- 
ure hy the energy method is developed, cased on substan- 
tially the same assumptions as those employed "by Wagner 
and Bleich. Problems treated in detail are: the stress 
and strain condition under St* Tenant twist and in twist' 
with axial constraint; the concept of shear center and 
the energy method for prohlems of elastio stability. 



1. CENTER 0? ROTATION AND WARPING UNDER 
ST. TENANT TWIST; SHEAR CENTER M, TWISTING 
WITH AXIAL CONSTRAINT 



Figure 1 shows the median line of an open section 
with several essential symbols. A reotangular system of 
coordinates x, y, z is passed through the oentroid S 
of one end cross section of the correlated column; z, y 
are arbitrary centroldal axes of the cross section, and z 
the column axis; axes x, y, z aro to form a right-hand 
system; i.e., the z axis in figure 1 points toward the 
observer. We arbitrarily fix a direction of rotation and 
thereby r.llocate to each point of the section center line 
a circumferential coordinate u. Suppose the direction of 
the positive sense of rotation indicates the positive tan- 
gential direction t; at right angles to it, the posi- 
tive normal direction n points toward the right, as ob- 
served when looking In t direction. The distances r t 

and r n of the tangent and of the normal, equal the dis-» 
tance r of the particular . point from the centrold S. 
A lino drawn from S in positive n direction indicates 
tho .positive rj. direction. This and the x direction 

form on angle a to be measured in ' the positive sense of 
rotation (rotating from +x toward ^y). If the t di- 
rection relative to has, say, a positive sense of ro- 
tation, r t is positive according to the preceding nota- 
tion. Correspondingly, r n is counted positive when n 
rotates positive in relation to S. 

Consider a column, as in figure 1, under the effect 
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of twisting moments I applied at Its tiro free end cross 
sections In conformity with St. Tenant*s theory of twiBt. 
The cross— sectional shape is preserved and all stresses 
and strains other than the angle of twist qp are unaf- 

fected "by z. The "reforred" angle of twist <p' = ^ is 

TP 

constant. With & = g (i+p) as 8near modulus and = 

j U- 8 3 as twisting strength"*"^ wheroby U is the devel- 
oped contour of the center line of the section and s is 
the wall thickness, it is T = G Jtp (p 1 . The linear dis- 
tribution of the shearing (twisting) stresses Tj over 

the wall thickness in any cross section is exactly the 
same as in a small rectangular strip. In particular, the 
shearing stress in the center line of the section is zero - 
as a result of which there is no angular change "between it 
and tho elements of the surface (longitudinal fibers). 
Every element of the modian line of the section remains 
perpendicular to the correlated fiber, which remains 
ptraight » according to the linearized theory. Then since 
the fibers, depending on their distance from the conter of 
rotation, are differently inclined, the cross section does 
not remain flat. If the rotation is, for instance, about 
an axis passing through S, the displacement of a point 
in tho cross-sectional plane is given with Y = r cp (fig. 
2); the component in tangent direction is accordingly: 

V t = r t <P (1) 

There being no angular change between the circumferential 
and the fiber element, the displacement is: 

Y = — -* + r- = 0 (2) 
dz du 

where W is the cross-sectional warping (positive in di- 
rection of the positive z axis). With the introduction 
of the unit warpin g w (of the dimension of an area, fig. 
1) as 

U=Ud 

1), 



For variable wall thickness Jj = -5* / a 3 du. 



■i-f 
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W = - qp' w (3) 



equation (2) (with equation l) gives - while observing ' 
that according to bb sumption , W, <p 1 , w are to he inde- 
pendent of z , = r t , 2 ^ that is: 



T 



du (4) 



'0 

If the rotation instead of about S is around any 
other point of rotation 33 with the coordinates and 

y D , through which a coordinate system x, y is placed 

parallel to x, y, it correspondingly gives 

V t = r t cp (5) 

W = - Cp 1 w (6) 

, u 

and w = w A + / r f du (7) 



J; 

ie f 

r t = r t + x B cos a + y 8 sin a (8) 



r o T / r t 
'o 

Figure 2 illustrates the following relations: 



and 



hence 



i.e. , 



or 



du cos a = dy, du sin a = - dx 

w - w 0 = w - w 0 + x B (y-y 0 ) - y s (x-x 0 ) 

w = w + x 8 y - y B x + K (9) 

w=w-x D y + y D x + K (9») 
__ . 

Strictly speaking, these relations are valid for the 

center line only; w = j r^ du + r n n is more exact. 

Since n, at the most, is oqual to s/2, the croBS-sec- 
tional warping r n n, superposing themselves on the clr- 

u 

cumferential warping / r t -du, can he neglected. 
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The unit warping w and w for one and the same 
section differ accordingly by one linear expression In z 
and 7. In "both cases the same stross condition prevails 
and the- reciprocal strains are completely the same. Tho 
difference in warping lios merely in tho relation to dif- 
forcnt reference pianos, which aro perpendicular to tho 
momentary axes of rotation of tho twisted column. Accord- 
ing to St. Vonant's theory of twist, tho position of tho 
axis of rotation is, on principle, undetermined bocause, 
since all floors remain straight, each fiber and each one 
parallel to it, may serve as axis of rotation. Fixing 
tho axis of rotation - say, through hinges - and prescrib- 
ing tho warping of any point, establishes the warping of 
all points if the load in the end sections is applied in 
accord with tho theory. But the axis of rotation can 
also "00 established, as seen from equation (9 1 ), "by pre- 
scribing the warping of throe (not lying on one straight 
lino) points arbitrarily, because threo of such equations 
suffice for the correct solution of tho unknown factors K, 
xt), and yjj. She froquontly entertained opinion that in 

St. Tonant's twist tho axis of rotation would always have 
to pas3 through ono definite point, tho center of shear M, 
i3 theroforo untenable. 

Since the concopt of shear center is being used re- 
peatodly in this report, a brief discussion of tho formu- 
las dofining its position is givon. Putting p. thin-walled, 
open-section column under tr?„nsvorsc load results, in gen- 
eral, in twisting in addition to warping. Warping is not 
accompanied by twist (definition of shear center) if the 
transverse force passes through tho center of shear M. 
Assuming linear distribution of tho bending stresses in 
trrirt-frco bonding and dofining the shoaring stresses due 
to transverse force in the usual manner from oquilibrium 
conditions, the statoment that in this case tho transverso 
force rolativo to the shoar center may have no momont, 
lead3 to the subsequent definition equations for the 
shear center : 2 ' 

it 

/ x w* dr = 0 and / y w* dP = 0 (10 ) 

'J p • 



'Tor the derivation of these formulas as well as of those 
used for integration of arbitrary sections, see the arti- 
cle entitled: "Shear Center of Thin-Walled Sections," by 
W. Luckor, oublished by the Static Test Branch of tho 
D.V.L. 
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Here x, 7 denote any coordinate - system passing through 
S , and the unl"t -warping*-, w*, ..referred to shear center 

M are given with "" 

.u 



■J V 

o 



w* = w 0 * + / r f .* du (11) 



whereby r t equals the distance of the circumferential 

tangent from the shear center (fig. 2). On the basis of 
the coordinate system x, y parallel to x, y placed 
through the arbitrary (to he chosen properly) reference 
point 0, figure 2 indicates 



Prom 



f ollows 



r t = r t * + x M oos a + y M sin a (12) 



u 

= i fi + / f * du (13) 



* _ 

r = w 



- x M y + y u x + K (14) 



The introduction of (14) in (lO)_then gives the two equa- 
tions (15) for tho coordinates xj( y^ of tho shear cen- 
ter M 



X M J x - Fit J X y= J J » dF 
- *M J xy + yil J y = -*J x w dF J 



(15) 



While no particular importanco attaches .to M in the 
oase of St. Tenant's twist, it plays a significant part 
in .twist with axial constraint, where cp 1 = dtp/dz is no 
longer constant. This case occurs, for example, if the 
column is clamped at one end or when, other than the end 
moments, individual moments or distributed moments. (m z ) 
are. applied in addition. Then the axis of rotation always 
passes through the shear center as shown later (equation 
18). All other fibers bond, because the warping V is 
not constant as a rosult of the changeability of qp ' and 
strains ** = /3a and stresses a? =.Ee*., occur in 
consequence in column dirootion. Inasmuch as the stresses 
<x* themselves then aro variable again, it simultaneously 
results, for roasons of equilibrium, in shearing stresses 
T* constantly distributed ovor the wall thickness, cours- 
ing in tho direction of tho column axis and of the circura- 
f eronco. 
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Unless the column length Is abnormally small, these 
shearing stresses T* remain small compared to the axial 
stresses a* , so that the related equations Y, in accord 
ance with the usual assumption of the, elementary "beam the- 
ory, may he ignored. The warping W of open sections is, 
therefore, under the assumption of preservation of cross- 
sectional shape, to be computed in the same manner as in 
pure St. Tenant's twist, where no slippage in the section 
line occurs either. From the partial differential equa- 
tion (2) circumscribed for T^.* and W in conjunction 

with the unit warping w defined by equation (11 ) fol- 
lows: 



hence 



W* = - qp" w* + f(z) 



<7*=Ec*=e|^ = - E(cp" w*+f«(z)) (17) 
o z 



As the twisted column is to transmit neither axial 
load nor bending moments, the axial stresses a* for each 
cross section must form an equilibrium group; that is, if 
F is the cross-soctional area 



J cr* dF = 0 , J x a* dF = 0 , J~* 7 cr* 



dF = 0 

'F w (18a, b,c) 

Since f'(z) is constant over F, and x and y are cen- 
troidal axes, the equations (10) follow direct from (l8b,c) 
which proves the coincidence of center of rotation and cen- 
ter of shear in this case. Equation (18a) gives: 



Compliance with 



f «(z) = cp" \ I w* dF (19) 
J 



J" 



w* dF = 0 (20) 

'F 

which is always possible with suitable choice of w 0 * in 

"(11) or of Z in (14), affords the following simple rela- 
tions : 

W* = - <P' w* (16 ') 

and a* = - E qp" w* (17) 
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The shearing stresses T* follow from the equilibrium 
condition f qj?= f orcos-in the , z „. direction at a small ele- 
ment s du dz: .... 

+ B |ll = o (21) 

a. u 3 z 

The shearing stresses t* d£ form a force couple of the 
magnitude - E 0 tp" 1 , whereby tho term 



0* = / w * 2 d] 



dP (22) 



is doslgnatod as warping strength relative to M. 4 ^ Tho 
twisting moment of the shearing strossos Tj being givon 
with G Jj qp' , tho term for the total twisting moment 
roads: 

T = 0 J T cp» - B 0* qp MI (23) 

Despite tho fact that the shearing stresses t* aro usual- 
ly substantially smaller than the shearing stresses Tip, 
both share approximately aliko on the total twisting moment, 
according to numerical calculations. This is due to the 
fact that tho lover arms of T* are of the order of magni- 
tude of tho section dimensions, while tho lever arm of Tip 
is oqual to 2/3 of tho wall thickness. 5 ) 

Tho calculation of tho angle of twist along the z 
axis for prodoternined individual twisting moments M z and 
distributed moments m 2 hinges on the integration of tho 
differential equation: 

- £2 = E 0* <p 17 - G J* «p" = m z (24) 

4 ^Since twisting with axial constraint is at times termed 
"flexurr.l torsion," the quantity G is also designated as 
toralon-ben^lag constant (G"bd or Cjjip (references 1 and 2). 
But in thq above case, the term "torsion bending" is mis- 
leading since twisting moments, but no bending moments, are 
transmitted (equations 18b, c). 

^This is roadily proved on a simple oxaraplo of tho I sec- 
tion. The axial constraint causes tho flanges to warp and 
act like individual bars; that is, have linear distribu- 
tion of a* and parabolic distribution of T* y while the 
web discloses no further stressos aside from Tip. 
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The limiting and transition conditions depend upon 
tho warping and stresses at the bordors; they are obtained 
from (15 1 ), (17'), and (23). With complete constraint, 
for example, we have cp ' = 0; with no constraint, it is 
qp» = 0, 

2. STRESS AND STRAIN CONDITIONS AT TWISTING FAILURE; 
ENERGY TERMS Aj AND A a ; DIFFERENTIAL 
EQUATIONS FOR THE STRAIN QUANTITIES £ R , T| e , «P 



To gain an insight into tho forces involvod at twist- 
ing failure, consider a contrally loaded open-section col- 
umn (fig. l) under the effect of the critical compressive 
force P. In nontwisted condition each fiber is under tho 
same compressive stroes; that is, in the critical case the 
twisting-failure stress is = P/F. While CT D and, con- 

sequently, the compression of tho column represent finito 
quantities, tho strossos and strains incurred on transition 
from tho straigh t to tho twistod equilibrium position, must 
be lookod upon as small quantities of tho first order of 
smallness. To them, the same relations, expressed in the 
preceding chapter with equations (I6')i (l? 1 ), (18a), and 
(19) to (23) apply, except that now the temporarily un- 
known center of rotation D is no longer coincident with 
tho shear centor M, since the statomont contained in (18b, 
18c) no longer holds true. The torque of the shearing 
forces in each cross section with respect to the axis of ro- 
tation D can be computed from (23^. where D substitutes 
for Li; accordingly, instead of C 

3 = / w a dF (25) 

whereby (?) and (9), respectively, are valid for the unit 
warping w referred to center of rotation D, and w 0 
and K nr-.st be so defined that 



/ w dF = 0 

For reasons of equilibrium, an anti-moment related to the 
external load is necessary; this is accomplished through 
the slope of the longitudinal fibers toward the axis of 
rotation at angle f qp' , so that the components of the 
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dompresslve forces <T D dff perpendicular to the ails of 
rotation, yield a rotating moment of 



ffj) dl r cp 1 f = cp 1 °"d Jp 

'TP 

about the column axis. The equilibrium of the momenta 
about the' axis of rotation then affords the differential 
equation: 

E 0 q? 7 + (0" n J p - & J T )qp" =0 (26 > 

from which with 

q> = a sin V 1 < 27 > 
follows the critical stress, as: 

& J„ + J o 

°b '—^ < 28 > 

J P 

But this formula does not as yet allow the calcula- 
tion of the critical stress Ojj, since the position of D 

remcine unknown. For the most general case of lack of 
cross-section symmetry, the_ two other conditions for mo- 
ment equilibrium about two cross-sectional axes are still 
needed. In twist with axial constraint, these equilibrium 
conditions lead exactly to the^shear center (equations 
(18b, c)); here, however, the 5 substituting for (7 must 
be so distributed that their moments on the assumedly iso- 
lated part, are in equilibrium with the moments of the ex- 
ternal forces due to the strain. Both moments are of the 
first order of smallness - the first on account of 5, 
the other on account of the lever arms. In Wagner's for- 
mula, whore D is equated to II, that ie,„tho same as 
(28) if D is replaced by M (that is, 0 by 0* and 

Jp by Jp ), these oquilibrium conditions are violated. 

Intending to resume these equilibrium conditions in 
section 5, we now proceed in a nimplo manner to the deri- 
vation of the differential equation of twisting failure by 
means of the onorgy method. To this the energy of the In- 
ternal and external forces (that is, tho total potential) 
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must "be expressed through the strains, after which the 
differential equations - the mathematical expression for 
consistency of stross and strain on the olement - thon 
follow purely formally in accord with the rules of cal- 
culus of variations. 

By A^ is meant tho difference between the strain 

energy of tho twisted and the straight equilibrium pat- 
terns under critical load P, and by A a , the energy of 
the critical compressive force P on transition from 
straight to twisted equilibrium position. Both expres- 
sions must - since this relates to a stability problem — 
be quadratic functions of tho strain quantities and, con- 
sequently, of the second order of smallness. Since the 
twisting process takes place of itself, i.e., without en- 
ergy input, the total potential of the internal and ex- 
ternal forces A^ - A a retains the valuo 0 of the ini- 
tial condition. Of course, the validity of Aj - A a = 0 
is contingent upon the correct strain in the onorgy terms. 
Por all other geometrically possible strains below the 
critical load, the energy required is A^ - A a > 0 , the 
very fact upon which the nonappearance of thoso strains is 
based. Por actually possible strains, i.e., compatible 
with tho oquillbrium conditions, A^ - A a is a minimum of 
value 0. This dual statement supplies the critical load - 
tho so-called stability limit, as well as tho strains (do- 
finod up to an arbitrary factor). Expressed in the form 
of the calculus of variations, the minimum requirement 
reads 6 (Aj - A a ) = 0. It always embodios the other dic- 
tum Aj - A a =0, as can bo proved. 

The strain energy A^ absorbed by the column of 
length I at twisting failure, is: 

z=l 

Ai = 1 J' {l € a dP + G J T qp' 8 } ds (29) 

2 = 0 



The assumptions we shall make are as thoso made for twist- 
ing with axial constraint. The cross section retains its 

shape, 6 ) so that tho strain may bo expressed by three quan- 
__ . 

'No floxural loads beine applied porpondicular to the sur- 
face of the column, the assumption of preservation of croos- 
sectional shape holds with great accuracy, failing only when 
tho longths or wall thicknossos are very small, as then tho 
walla nay bucklo. 
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ties: For example, the angle of twist cp and the dis- 
placements | 8 and T\ g - of the - centroid S. relative., to 

tho position of the center of rotation D, . which are to 

ho measured in the direction of the x and 7 axes of 

the cross sections of the untwisted column. Then. 

W = - cp' w (30) 
7. jf.-,-* (31) 

replace (16»), (17 1 ), and (20) 
and 



J. 



i 

w dJ = 0 (32) 

F 

whereby (32) voices the equilibrium condition essential to 
tho stability problem, namely, that the axial force (-P) 
remains unchanged. Equation (9) herewith becomes: 



- 1 !\ 



w = w +x_ y - y_ x - A / w dF (33) 
with I" 

Is = " 7 B "P. Tl s = + x 8 cp (34) 

then follows from (31) the important relation for the 
strain distribution: 

/■■» 

* = " Is" ^ - n B " y - cp" (w - \£ w dP) (35) 

after which a simple computation and the introduction of 
(35) in (29) gives the strain energy Jl± : 

z=l 

its 



-»/< 



B J y | S " B + 2E J xy | 8 " n 8 " + X J x T1 B T 
a=o . + 2B Hy Is" cp" + 2B H x Tl s " «p" + 

+ E 0 qp" 8 + d J T cp ,a }dz (36) 
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Herein: 

R y = / x w dP (37) 
R x = / y w dP (38) 

y w a di- - j (y w dp) (39) 



c = 



Che quantities R y , R x , and 0 are dependent on the 

ahapo of the cross section only. In conformity with the 
concept of stat ic moment R x (and Hy) are termed "curv- 
i ng moment for the x. axis" (for tho y axis, respective- 
ly), C the "cu rving strength relative to S" (cf. C*, 

equation 22). For symmetrical sections, for instance, the 
course of the unit warping w is ant i symmetrical , if the 
point of the section center line which meets the axis of 
symmetry is coordinated to zero warping (equation 4). 
Thus, Ry = 0 for symmo try with the x axis, and R x = 0 
for symmetry with tho y axis. As z, y are centroidal 
axes, these relations are applicable to any other determi- 
nation of tho arbitrary constant w 0 as well. 

Tho energy A a of the external forces - that is, of 
the compressive force P - due to contraction of chord on 
Duckling of the fibers, is computed as for the "buckling 
Euler column, since as a result of the rotation about the 
axis of rotation D all originally straight fibers change 
to "Piano curves. With | and T| as displacement compo- 
nents of any cross section point, it is: 

z=l 

A a = a D ! / / U ,s + Tl' a ) dP dz (40) 
z'do J P 

According to Suler's theorem of kinematics, the motion of 
a point is built up from tho motion of any reference point 
and the rotary motion around that point. Por S as ref- 
erence point, it thon is. 



t = 5s - y <P. Tl = Tls + i <P 



(41) 
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Consequently, since z, 7 are oentroidal axes, 

z=l "" " " '■ ■ 

= ?L i J t {iT(Ib ,B + T1s ,B ) + < J x + Jj)*' - } d « < 42 > 
z=o 

According to tho rules of calculus of variations: 

8 (Ai - A a ) = 0 (43) 

then affords in simple manner the three differential equa- 
tions for tho strain quantities £ a , T\ g , cp: 



[l! + S>^ff a IV + 3E>y^ 7 + O p P.e a " = 0 " (44a) 

E^ y '4 IT + B J T Tl. a IV + S^g^ 17 + Op F tl a » = 0 " (44-0 ) 

I ^<t fl 17 + E^%, 17 + E 0 <p 17 + (Op Jp - G J T )q>" =0 ^(44c) 

They state that the forces following from (30) must he in 
equilibrium at each column element (six equilibrium con- 
ditions in similar manner as this is expressed for a Euler 
column "buckling in the xz plane hy the differential equa- 
tion E J y | 17 + P |" =0, which, with observance of the 
strain law B y = - E J y J" follows from the equilibrium 

conditions ^— 2 = Q, and = P | . 

oz oz 

Because equations (44) are coupled together, the three 
quantities | fl , T) B , and cp usually occur concurrently; hut 

if symmotry prevails, it is otherwise. With symoB try to 
tho x axis, for instanco, equation (44a) is independent 
of (44b ) and (44c), because J xy = 0 and E y = 0. In this 

case a deflection g 8 , independent of T) B and cp Is pos- 
s i bl e • 

In tho simples.t case of double symmetry, tho throe 
equations and with them tho strain quantities £ 8 T| B , and 
cp thomselvos, aro lndopendont of each othor. Accordingly, 
thoro are then, three different buckling processes, of 
which, of course, only the one with the lowest critical 
compressive stress is of practical significance . The crit- 
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ical compressive stresses 0" Dl = ff yf <jj ig = ff x , = cr s 

i g = a! sin 3L8 f r\ Q = ajJ gln 113, qp = a 3 sin ^ (45a, b,c) 

at 

CT 7 = nr — Z . *x = -3" — " (46), (47) 

IE IE 

and 

G J T + ^ E 0 

ff — 1 (48) 

8 J P 

Under tho critical twisting-failure stross, according to 
(48) tho column thus twists about an axis of rotation, 
passing through tho controid S, Eor an I section, for 
example, it means, with g = 2.6, if h is tho hoight of 

tho web, and b is tho width of tho flange: 

ff hj b 3 + Q. B l g (h ± 2b? S a \* (49) 

l B 4b 3 + 12 b h 3 + 2h 3 

To make twisting failure possible practically, tho I soc- 
tion must not only havo the proper thin wall, but its flange 
width must oxcood tho hoight of tho wob; for tho condition 
of smaller a B then tho Eulcr stross for buckling about 

tho wob axis, roads: 



(£) + 0.468 (h + 2b) s a V £ 1 (50) 
W b 6 



Eor sections without axial symmotry but disclosing 
point symmotry (for example, Z sections with equal legs) 
(oquations 44; also aro not coupled because H x = R y = 0 , 

rogardloss of whothor tho controidal axes x,y ar.o princi- 
pal axes of inertia or not. And tho twisting failuro again 
occurs about an axis of rotation through S undor criti-* 
cal compressive stross <T a , according to oquation (48). 

In tho most gonoral caso of asymmotrical open sec- 
tions, none of tho stresses cited in equations (46) to (48) 
is a critical -compressive stross. The critical values 
CTjj x . CTjj a t (Tj) 3 then follow as roots of a cubic oquation. 
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With simple symmetry thoir place is taken by a linoar and 
a quadratic equation.. 

Bloich's article lacks the essential term o~p(J x +Jy) 
<p" in the expression for A & , though the expression for 
the strain energy A^ is as to contents in agreement with 
(36 X Closed formulas for the quantities Hy, E x , C of 
the typo of (37), (38), and (39) are not given; tho con- 
copt df unit warping itself is not employed. Their inves- 
tigation is restricted to contoured columns built up from' 
■flat plates. They proceed from the curvatures of the in- 
dividual plates expressed with throe strain quantities 
l S f T| B , CP, and to which the conventional theory of beam 
flexure is applied. Tho strains (?) are mado proportion- 
al to those curvatures and a linear distribution is postu- 
lated along the straight pieces , of which tho cross sec- 
tion consists. With n plates it affords for tho calcu- 
lation of the n unknown strains at plate center, n - 1 
transition conditions and one equilibrium condition (aqua- 
tion (32)). The calculation of those n unknown strains 
must bo carried through first for oaoh section. And even 
then the expressions for H y , H x , and C (given for some 
sections by 31eich) follow only after tho dovelopmont of 
tho expression for A t . Tho lack of torm 0"i)(J x +Jy) <P 18 
in tho cxpro3sion for A a rosults in tho completo absonce 
of Oft in Bloich's equation, which corresponds to equa- 
tion (44c). 

3. SPECIAL CASE 03? SECTION SYMMETRICAL WITH 
X AXIS, CHANNEL SECTION AS EXAMPLE 

For symmetry with the x axis, it is: J xy = 0 and 
R y = 0. T^ls releases (44a) from (44b) and (44o), which 
are coupled through R x £ 0. The solution is again ef- 
fected with (45) but now ' a 8 and a 3 are no longer inde- 
pendent of each other. Equation (44a) becomes: 



e jy e 8 IV + ct d j e 8 " = 0 



(51) 



with which (45a) becomes: 



(J. 



= 0", 



■n 



8 E Jy 



(52) 



y 




Erom equations (44b, c), now written: 
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H J x T1 S IV + ! Bj cp 17 + ff D E T] B " = 0 (53a) 
E IL,. Tl a 17 + E 0 qp 17 + (O'd J p - G J T ) <p" = 0 (53b) 



follow tiro homogeneous linear equations for the unknown 



(54a) 



a a an d 

[jjSHxjaa + ■ 0 + G J T - ff B Jp J a 3=0 (54b) 

which have a solution other than zero only when their de- 
terminant disappears: 



TT 



E H, 



= 0 



(55) 



This condition affords a quadratic equation for the 
two critical compressive stresses Og a and °i} 3 ' Abbrevi- 
ating (47) and (48) with the added abbreviation (dimension 
of a stress ) 



w 3 E E x 

3C = ~"S 



1 /Tj! 



(56) 



equation (55) can be written as: 



= cy D a - ct d (a a +CT x ) + cr a cr x - p x a = 0 



that is: 



(57) 



(58) 



This relation can be represented by Mohr'e circle 
(fig. 3) and thus enable a simple graphical solution of 
Oq from CT g , CT X , and p x . The smaller root is seen to be 
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consistently smaller, the greater root consistently great- 
er than cr_ and o"_ . The strains correlated to c-n and 

ctjj are given with equations (45b t c), from which follows: 

The cross sections accordingly twist about the centers of 
rotation Da and D 3 lying on the axis of symmetry and 
whose coordinates x Da and xjj 3 , after insertion of <j- Da 

and (Tjj 3 in one of (54), give: 

*D = - I s = -Jlllh /5 __ P *_ (60) 

If 

or else g 

0 J T + itg E C - ffp J x 
x = - fJL L 1 = / I£ rJL.Z_T£ (51) 

Por p_ ^ 0 it Is c p ^ cj^ and <Tjj ^ ^s* the value 
for xp is, in consequence, always finite. In other 
Tvords, re obtain the important rosult that on a symmetri- 
cal open section a buckling in direction (y) perpendicu- 
lar to the axis of symmetry (x) is always accompanied by 
column- twist. For very thick-walled sections, of course, 
cr 8 is very great on account of the groat twisting strength 

so that o~2 a » o* z and o~d 3 » o~ fl . This is recognized from 
Mohr's circle or else from tho approximate formulas valid 
for groat ratio O'b/'Jx 5 

P 8 




a^-a x -5j- (62) 



3 s 



following from (57) or (58). In this case, for au a * 
a z the centor of rotation Dg is extromoly remote from 
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tho centroid, while for very high twisting-failure stress 
a$ 3 « or B , the center of rotation D 3 is almost coinci- 
dent with the centroid. On sufficiently thick-walled, sym- 
metrical open sections, the "buckling in direction perpen- 
dicular to the axis of symmetry is therefore practically 
free from twist. 

Bleich, lacking tho term -ojj Jp in the determinant 
(55), reached instead of (57), the linear equation: 

<*x - bleich) *s - Px° - 0 (64) 
7) 



with the single solution 

CT Bleich = ^x ~ cf 



Px 8 



(55) 



s 

This formula which, with observance of 

°x " Jx/Jp " *P 

(oquations 47, 56, and 15) nay also be written as 

s 



(66) 



agrocs with tho approximate formula (62), valid for vory 
small ratio a x /cJ B , and supplies for that reason useful 

values for the conventional sections and lengths used in 
structural engineering. But for the thin-walled sections 
customary in airplane design, the discrepancies between 
°Bleich an ^ °D ma y become quite considerable. Prom the 
theoretical point of view, it is even more essontial that, 
according to Bleich 's formula for A a , one of the criti- 
cal compressive stresses always be lost. According to that 
theory, for instance, twisting failures of point symmet- 
rical sections should be impossible. 
_ 

/ Bloich , s equation of the critical load, written with our 
notation, reads as follows: 

H_ a \ it 8 I J 
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The "boundary conditions for the case corresponding to 
(58) are, as subsequently established from (45) for both 
cases (z = 0, z = I) 2 T1 B '«= 0, TV B n «* 0 and q> = 0, 
cp =0. It means, according to (31), that 5*0; that 
is, complete axial nonconstraint obtains. For it, W = 
- cp 1 w with <p' jfe 0, according to equation (30). The 
exact realization of these boundary conditions, even on the 
simplest cross-sectional forms, would require special sup- 
ports with several knife edges, in order to enable the in- 
dividual parts of the support to follow the cross-sectional 
warping. 

TTith complete axial constraint (say, by welding to 
rigid support plates, etc.), the following formulas replac- 
ing (45) are appropriate: 

^a^as^l-cos "f—^O" c P=a 3 (l-cos (68) 

At both boundaries wo havo cp = 0 and cp 1 = 0; that 
is, TT = 0 , according to equation (30), and in addition, 
"H B = 0 and T\ e * = 0. As cp =j£ 0, 5 likewise is of ^ 0, 
according to equation (31). Everything else follows in 
tho samo mannor as beforo; in equations (54) to (58) l/2 
replaces l . 

For other boundary or support conditions, it is recom- 
monded to obtain by tests a substitute length which re— 
pl.icos I in (58). Tho substituto length lies always bo- 
twocn 1/2 and I. 

Tor illustration, we repeat the formulas for a channel 
section. Tho solution of the unit warping, of warping mo- 
ment IL^, and warping strength 0, is carried out In soc- 
tion 8. TTith a as web height, and b as flange length, 
it is: 

CT - IlfLJ! a 8 (a+6b) a _ it 8 g b 3 (2 o+b) ' ( . g) (7Q) 

°= - I 8 12(a + 2b)' ~ IT 3(a + 2b) a (69) ' ^ 

■ _ a a b 3 (2a a +15ab+26b a ) + ±- |(a*2b) 3 a" I 8 , % 

Q- _ n_s : tt ■ e . (71) 

8 1 (a+2b) |a s (a+6b) (a+2b) + 4b 3 (2a+h)J. 

p x = j£= g a^ b 8 (a+3b) (72) 

1 (a+2b) y5~^ a (a+6b) (a+2b) + 4b 3 (2a+b)| 
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In the special case of a = b,. illustrated in figure 4, 
we obtain (with | = 2.6): 

a _ = J- Tf3E a = 0.1944 — — (73) 
X 36 (l/h) 8 (l/h) 8 

a. = [ 0 .4345 ♦ 0.04251 (|f J^ljr (75) 

Pl = - 0.2680 ■ (76) 

Then .ojj follows, according to equation (58) or with the 
aid of Ifohr's circle. We have included in figure 4 the 
curves for Wagner's fornula 8 ): 

G J T + =5 E. C* 

(77) 

J P 

as well as Bleich's formula (65) or (67) for comparison. 
For a = b, is obtained: 

9 U = 0. 02239 + 0.01466 (*) (£) \ (*B) 

and 

■ta.i.h-0-"** [l O^MM—I jtt^ ,„) 

0.4345+0.04251 (f) (J) J (l / 1>) 

The discrepancies between a^j and cr^ as explained 
in section 6, increase with wall thickness and length, 
while tho c" if f erence e between °Bleich an ^ °D are S T0& ^" 

8 ^C* is identical with Wagner's C^^. Ho then addod (dis- 
regarded hero) the term &bdn connected with the transverse 
warping, but which servos no usoful purposo unloss it pro** 
cisoly pertains to vory short columns of unflanged anglo 
section, T soction, or cruciform section. 
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oat for small length and wall thickness. The curve of the 
twisting failure stresses ct d cannot approach the Eulor 
curve Td^) "for buckling perpendicular to the axis of 

symmetry, except asymptotically, according to equation (62), 
although it can intersect the Eulor curve (a y ) for "buck- 
ling in direction of the axis of symmetry. This is the 
caso for so much shorter lengths as the section wall is 
thicker. 

4. THE GENERAL CASE 



If the cross section reveals neither axial nor point 
symmetry, the warping moments | 8 , T| B . cp do not disappear 

and equations (44a, h, a) are, as a result, couplod. In 
equation (45) 

t TTZ ^ TTZ TTZ 

E s = ai sin -y , Tl 8 = as sin y, cp = a 3 sin — 

therefore the unknown coefficients a lf a a , a 3 , are mutu- 
ally dopendent and result in throo linear, homogeneous 
equations, whoso determinant must disappear for the criti- 
cal caso. This relation supplies a cubic equation for the 
three potential trri sting-failure stresses CTj) , 0"D a , 0j) 3 . 

of which, of course, only the smallest retains our intoroBt 
The insortion of the ff D values gives three different ra- 
tios ai/a 3 and as/a 3 oach, which define the position 
of the three centors of rotation Di , Da, and D 3 , as is 
soon from the following equations: 



U - aT * = * 5% * = + y D 80 that 7 D " aT 



Tig = — cp = + x B cp = - xjj cp, so that x D 



a 3 J 



(80) 



With the ahhreviations (46), (47), (48) , and (56), as well 
as 

7 1 J* Jp 



*xy = Jf C82} 
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the disappearing determinant 



TT 

7 



TT 



35 J, 



It e ^ 



fa ■ H x 



2| E C + G J T 



ff D J p 



= 0 



may "be written in the form: 



P„ 



'xy 



P: 

cr x - crjj 



CT s - CT D 



= 0 



(83) 



Equation (83) ia the well-known "secular equation" which 
reads : 

a D 3 - ^D B (^ 8 +o x + ^y) +<Ttj(<t s <7 i + Cfy + cr x or y . 



- ?xy 8 - Px 3 - Py 3) - (°b¥ 7 + 2 PxyPxPy - 



- P a 



a - p 8 a - p 8 a ) = 0 
p y - r x y 



(84) 



All roots are real (symmetrical determinant 1 .) That 
they are always positive as well, is seen from the follow- 
ing. The integrands of A^ and A n in the strain quan- 
tities £ B , Ti Q , cp and their derivatives, respectively, 
"being quadratic and homogeneous, Ai and A a "become 
through (45) homogeneous quadratic functions of the coef- 
ficients ai, a 2 , and a 3 . The result is: 

±i - I £ [{4.* J y a! 8 + 2 H| E j. 



"xy a i 



+ 2 



TT 3 



a a » 

^ > *y a 1&3 + + ((J J T + J| EC) a 3 8 J 

= $ (a x , a 8 , a 3 ) (85) 
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and 

*a » V-*£-V +•*■■■). + Jp^a] = 

= CT D +*(a Xi a a , a 3 ) (86) 

Since tho difference A x - A a and with it 9 - Oj) i|/ 
shall "be a minimum, the three equations 

- „ * — = 0, 1=1,2, and 3 (87) 
o &i 

must do complied with. Obviously, these are identical 
with those obtained from (44) through (45), whose disap- 
pearing determinant leads to equation (84) and so, to tho 
throo values oj) , crn , op, and to the three ratios ai : 

aa * d ji To <b - oj) \J/ as homogeneous quadratic function 
of the a 4 , tho Eulcr formula: 

8(9 - a. MO 3(* - or,. 

±! ai + a 8 

d a x 8 a a 

3(9 .. 0"_ \|n 

+ — =—2 a 3 = 2(9 - Ojj *) (88) 

o a 3 

is applicablo; that is, the difference "becauso of tho 
throo oquationB (87) is 

9 - CT D \|/ = 0 (89) 

which in tho present case proves the correctness of the 
double statement Ai - A a = min = 0. Such a proof can he 

quite generally adduced when introducing 

Is = »i MO, T1 B = aa f a (a). «P = a 3 f 3 (z) (90) 

in place of (45), because 9 and \|/ are again always ho- 
mogeneous quadratic functions of a x , a B , and a 3 « The so- 
lution of equation (89) with respect to Oj) gives: 



9(a lt a a , a,) 



9 



\a, a,/ 



a D = IlZllJi*' "=» = ™* (91) 
v|/(a lt a a , a 3 ) ^ /ai os\ 

via' a s / 
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a i &q 

where, of course, the values — and — should he in- 

&3 as 

sorted (hut which follow from (87), not hef ore ) . Having 
thus defined \|/ and $ as positive, as 1b readily seen 
for \|/ from (86) and for v|r from its. physical signifi- 
cance, it follows from equation (91 ) that the throo roots 
Oft , CTtj , and OV, are all positive (i.e., compressive 

strcssos). 

Tlio premise A^ - A a = min = 0 can, moreover, "be re- 
placed "by the other demand: 

^ = ^L^rv^)J C92) 

therein | , T) , cp would "be the solution functions (un- 

knorn as yot) of the prohlem. The differential equations 
for thcso functions r.nd also of a D itself, are ohtainod 
from equation (92). The variation carried through, gives:' 

>|/fi$-fi)6\|/=0 (93) 

6 $ - | 6 \J/ = 6 ($ - \|/) = 8(Ai - A a ) = 0 (94) 

which proves the previous argument. 

The practical solution of (84) is "best effected with 
the graph for the cubic equation given as figure 5: 

uu 3 - 3u) e + p cu - q = 0 (95) 

Hereby: 

p = g s °x + ^s g y + ° x Py ^- Pxy a - Px a - P y a ^ g6 j 

°m 

= J a °± g y + g Pxy Px Py - Px j* g s - Py e g x - Px 3 g y (g?) 

0- m = — Z (98) 

First determine: <J g , <J X , Oy , then cr m , p, q, and read 
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the particular minimum U) value from tho chart. The 
twisting-failure stress then follows as 

ff D = «" CT m (99) 

Por very thick-trailed sections, the conditions resem- 
ble those of the special case discussed previously. One 
of the three roots, say, 0|q 3 , is approximately equal to 
the very high stress a a , while the other two roots fol" 
low from equation (84) divided "by o" s . The "buckling here 

is practically froe from twist. With x and y as prin- 
cipal a=cs, for example, (P TT = 0), it is oj), * <7 X and 

(J_ ~ (T 

d 8 . - V 

5. REPERENCE TO ANT EEPEEENCE POINT 0; TWISTING PA I LURE 
POH POECED AXIS OP E0TATI01T A 



The cross-sectional shauo Doing preserved, tho de- 
scription of the strains can be made with throe (on z do- 
pondont) sots of data. This time, howovor, the choson 
strain quantities include, asido from the angle of twist, 
tho di aplacoment s | 0 , Tl 0 of an arbitrary roferonco point 

0. 0 is to have the coordinates x 0 , y 0 and from the 
origin of a coordinato system x, y parallel to x, y. 
The unit warping w roforred to tho roal center of rota- 
tion D is oxprossed, according to figure 2 (of. equa- 
tion 9), as: 

. w = w + x 0 y-y 0 x+K (100 ) 

Hereby the unit warping W roferrod to 0 is given 
throush 

w = w 0 + / r" t du (101) 
•J 

Purthor, it is " 0 

&o = ~ ?o "P. Tl 0 = + x 0 q) (102) 

that is: 

T = - lo H x - Tl 0 " 7 - <p" w - J / w dP (103) 

i 
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The result is: 
z=l 

±i = i f ( E J 7 + 2E Jxy 5 0 " Ho" + 1 J x Ho" 3 



z=o 



+ 2E Ry £ Q " qp" + 2E E^. 71 0 " + 

+ E "S cp" 8 + G J* cp |3 Jdz (104)* 
S y = y x w dE (105) 

R x = / y t: IP (106) 

«/ 
E 

0 = / w a dE - \ ( u IP) (107) 
./ F V ' 

iloroovcr: ' ^ 

t = ! 0 - 7 «P. Tl = T1 0 + x cp (108) 

and 

?b = - y 0 - *s = " x o (109 ^ 

z=o 

+ 2E 7o | 0 « qp' - 2F x Q T) Q " qp ' }dz (110) 

Herewith follow from 6 (Aj - A a ) = 0, tho throo difforon- 
tial equations: 

EJ y e 0 1V + EJ iy T1 Q IT + E E y 0,17 + d D F e 0 " + a D F y 0 cp" = 0 

(111a) 



with 
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EJ xy ?o 17 + EJ x ^o 17 + El x ^P 17 + O33 F Tl 0 " - F x 0 q>" = 0 

-(HIT?) 

Ely £ 0 IT + Mt,. T) 0 IT + 30 <P IT + ff D 7 y 0 l 0 " - (T D F x 0 Tl 0 » + 

+ J p - 0 J T ) <p" = 0 (111c) 

These equations (ill) are more complicated than equa- 
tion (44), hut they can he considerahly simplified when 
introducing the actual center of rotation D as arhitrary 
rpferenco point O p because 

l D = 0, n D - 0 (112) 

They then "become: 

3 E y cp IV + a D F y D qp" = 0 (113a) 

3 Ej. <p IT - ff; P i 3 q>" = 0 (H3h) 

E C q> IV + (a D J p - fl Jrn) <p" = 0 (ll3c) 

whereby tho warping momontB: 



E, 



x 7 

and F 



= y z w dF (114) 



£- = / y5dl (115) 



F 

as woll as tho warping strength 0 depend on the (yet un- 
known) position of tho center of rotation D. liquation 
(113c) in in accord with equation (26), ohtained from the 
condition for momont oquilihrium ahout the axis of rota- 
tion of a column oloment. The equilibrium condition £Z = 
0 for tho forces in column direction is alroady mot hy 
equation (32) 

7 dF = 0 



J 



F 

thus oliminating tho second term in oquation (107) and 
leaving for 5 tho simplor oquation (25): 
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-a 



0 = j w a dP 



Equations (llSa.b) are tfce previously cited conditions for 
equilibrium of tho moments about the y and x axes in 
connection with the equilibrium conditions EI = 0 and 
ST = 0 for any column element. This is readily seen from 
the expressions for the bending moments and transverse 
forces. TJith equation (3l) the bending moments are: 

x a dE a - B Ey op" (116) 

and r 

B x = i y dP = - E E x cp" (117) 
and the transverse forces: 

U=Ut, 

/ * / • 

Qac = - j T s du sin a = + J f s dx = - E Ey <P»" (118) 

u=u D ' I 

and 0 

u=u p ' 

/' /" 
Q y = + / f s du cos a ^ + / f s dy = - E S x cp m (149^ 

u=u- " E 

when applying the relation: 



u 

T s = + E cp 111 J if dP (120) 
u=u. 



following from equation (21) and effecting a partial inte- 
gration. 

Writing the expression (9) or (9 1 ) into (114), (115), 
and (25), givos: 



7 = S y + 5 0 J xy . - y g J y = E y - x D J xy + y D J y (121) 



E 

*z ■ H x + *s J x - y 8 J xy = *x " *D J x + * D J xy < 122 > 
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and 0 = 0 + x D a J x - 2x D y D J xy + y D a J y 

- 2x D E x + 2y D H y (123) 
Then oquation (113) can bo eolrod with 



cp = a Bin ™ (124) 



yielding at first: 



^ = 4^, ^ = -4^ (125), (126) 



and 

0 J T + E C 

I 

D = = 

J P 



(127) 



Ey oliminating Oft from t»aeh of two equations , the 
unknown distances x^ and yjj of the center of rotation 
can he determined, after which insertion in one of the 
three (125) to (127) gives oft. Compared with (95) to (99), 
this method is extremely tiresome, since the equations for 
Xjj and y D are couplod and of the third degree. 

Admittedly, considerable simplification obtains if 
symmetry prevails. With Ry = 0 and y D = 0, for ex- 
ample, equations (126) and (127) become: 



a„ = - Hi Hz x * J *> (128) 

l a I x D 

and 

. &-J T + H| B(0 - 2x D R x + x D a J x ) 
0- = 1 (129) 

J P + * V 

Prom this follows as distance xj of the axis of rotation, 
from the center of gravity, the quadratic equation 

x . _ x a_5 »_ E = o ("o) 

K X 
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Owing to the term (~tt) (J^ Jrp I 8 , x D is not only de- 
pendent on the cross-sectional parameters given through 
the section center line', hut also on the wall thickness 

and the column length (parameter , where b 

denotes any cross-sectional dimension. 9 ) ) 

Effecting any axis of rotation A intersecting the 
axis of cymiaetry - by hingelike guidances along the whole 
"bar (fig. 6a), for instance - affords the following equa- 
tion: 

G J T + 1(0 . 2x A E x + x A a J x ) 

<r A = 1 (131) 

J p + I x A a 



which springs from equation (129), if x D is replaced "by 
x A , the distance of tho axis of rotation A from the col- 
umn axis passing through S. For, as is readily seen, the 
reaction forces appearing on the guides have no offset on 
tho equilibrium of tho moments about the axis of rotation 
A, and oquations (26), (113c), and (129) retain thoir va- 
lidity, if A is used instead of D. But now equations 
(ll3a,b) ere no longer applicable, because tho roactions 
modify the transverse forces in the individual cross sec- 
tions. Though the onor-y expressions Aj and A a (equa- 
tions (104) and (110)) remain the same, only <p may be va- 
ried in 6 (A^ - A a ) = 0, as cp now constitutes tho solo 
indopondont quantity. 

Figure 7 illustrates tho relation of critical com- 
pressive otross OA to distance x^ for a certain chan- 
nel section. The extreme points of the curve = f(x^) 

are at the sane time the zero points of equation (150), as 
__ 

'While reading the proofs of this article, a report touch- 
ing on the same subject by 35. E. Lundquist, entitled: "On 
the Strength of Columns that Fail .by Twisting," appeared 
in The Journal of the Aeronautical Scioncos, vol. IV, no. 
6, April 1937. Hef erring to H. Wagner (reference 4), 
Lundquist proceeds from equation (127) and postulates for 
the axis of rotation such a position as will cause o-q 
to be a minimum. This method leads to the solution of an 
equation system - (125) to (127) - (or of (128) and (129) 
for symmetry). 
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can "be proved from a differentiation. The two roots of 
equation (57) -therefore agree- with, the extreme value b of 
equation (131). This is directly connected with "the pre- 
vious formulation of the energy method in equation' (92), 
"because the numerator or denominator, respectively, of 
equation (131), is proportional to the respective internal 
or external work in establishing the sine formula (124). 

For x^ = Co, ct a is equal to the Buler stress <j x . 

This case can he enforced, for example, with guidances of 
the type of figure 6f>. The result of the force couples 
due to the guides is that the transverse forces applied at 
the -cross sections always assume the position demanded "by 
the equilibrium. On the channel section of figure 6b, for 
example, the transverse force Qy must pass through S on 

acoount of the central compressive force, which is impos- 
sible- without guides. 

Tho cases cited in figures 6 and 7 are of significance 
insofar as columns of open section are frequently used as 
stiffeners of the skin of stressed-skin structuros. The 
axis- of rotation is in this case subject to considerable . 
displacement toward th3 sheet which in addition furnishes 
an elastic support against twisting, resulting in a marked 
rise of critical compressive stress compared to the values 
computod hero. A nathonat leal analysis of the two effects 
on combined action of metai skin and contoured columns is, 
of course, quite difficult." , 



6. RELATION OP TWISTING- FAILURE STRESS Ojj TO 
WAGNER'S CRITICAL STRESS Ojj 

With center of shear li chosen as reference point 
in tho sense of the preceding chapter, we obtain, in place 
of equation (ill), three corresponding differential equa- 
tions for the displacements of the shear center T^, 
and tho angle of twist op. With the unit warping w* as 
deflnod in equation (11 ) and referred to' the shear center, 
C now replaces t», that is: 



= J w* a dF - £ ( / w* dF ) 
'" F 'F ■ 



(132) 
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Onlng^to equation (10), the momente Hy* and • E x * replac- 
ing Ry and E z disappear, which affords a certain sim-- 
plificatlon over equation (ill): 



dF = 0, 



y w+ dF = 0 



(133) 



Only tho symmetrical case is treated hereafter. For 
symmetry with the x axis, we ohtain with J xy . = 0, 

yjj = 0, the following equations from equations (111b, c): 10 ^ 
E J x T) M IT + 0- D F - Ojj F x M <p" = 0 (134a) 

E C* <?> IT + (ct d J p * - a J T )cp" - O-jj F x M T1 M » = 0 (134b) 

Using tho abbreviation 0"jj for Wagner's critical stress, 
according to equation (77), in conjunction with equation 
(45) for the critical compressive stress Op, we have: 



Xu 

°x — °E °D 7 ' 
1 P 



0\r - OV 



/ Xif a \ 



<0 + ^M a x = 0 



that is: 



l CT v +C7 M + / (^t-O'v) 8 + 4CJ T °}l -^-»T 

1 P 



(135) 



(136) 



or else: 



JL = I (A. + _1_ ± /Zl _ JLV 
ct d 2 la M <r x J \<j u <t x J 



°H ff x l p 



(136 ») 



a z are critical 



If x ^ 0, neither ajj nor 

stresses, as seen from equations (135) to (136 *)• The 
relations! 

10T 



Cf. equations (5Sa v b). 
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-D 



and 



following from equation (136 1 ) are illustrated in figure 8. 
It immediately affords the percentage of error incurred 
when tho critical compressive stress <7q 1b computed "by 

Wagner's formula (77) or according to Euler's (47). The 
error increases with the ratio, respectively , of O-^Ja^ 
or" o^/ojj aad with tho ratio ^14/1^*1 dependent ' on the 

section form as, for example , is shown for channel sec- 
tions against- aspect ratio a/h in figure 9. The dis- 
crepancies "between 0"d ani °"M' and oetweon and ct Xi 
respectively, are least when the ratio 



°X T 1 a 

* J P *X 



(139) 



is a minimum or maximum, rospoctivoly ; i.e., in tho case 
of vanishing wall thickness or column length, re spectivoly , 
for very groat wall thicknossos or column lengths. In 
theco casos tho approximate formulas: 



or 



% " 1 + £ (141) 

are applicable. The latter is in agreement with the ap- 
proximate formula (62) and with equation (67) for CT Bloich' 
For tho case Ojj = O" , 

is exactly valid. 







X 


i * a 






% 


i * a 

1 P 



(14Q) 
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If limited to values "below 1 for qjj/a" x (shorter col- 
umn lengths and vail thicknesses), respectively, and for 
°x/ CT M (greater column lengths and wall thicknesses), the 

error for ctjj conformable to (77) and (47), respectively, 
cannot exceed 100 percent, according to equation (142). 
The maximum error in the most unfavorable case is 97 per- 
cent for channel sections, according to figure 9, and 82 

percent for the case ^ = 1^, illustrated in figure 4. 

In reality the maximum error hero does not exceed 44 per- 
cent (fig. 4), "because the Euler stress or y for "buckling 

in symmetry direction becomes less than <7 M for fairly 
short lengths, bo that starting with such lengths the crit- 
ical stress would then no longer be computed by Wagner f s 
formula (77) but by the Euler formula (46). Such limita- 
tion of the maximum possible error exists, for example, on 
all channel sections with a/b > 0.73 (a_ < cr^). Figure 

8 shows the maximum values of GuZ^D ^ or channel sections 
with a/b > 0.73 as a dashed curve. Tor the rest, tho 
use of figure 8 in conjunction with Wagner's formula (77), 

is of real advantage only if tables for 0 , J , and xy 
arc availablo; otherwise, equation (58) gives quicker re- 
sults. 

■ * « 

It also will be noted that the Wagner stress 0"jj is 

directly obtainable from the enorey method through the 
arbitrary (Ritz's) formula: 

In = 0, T1 M = 0, cp = a 3 sin ^ (143) 

which, written in (104) and (110) for Aj and A a - cen- 
ter of shear M replacing reference point 0 - givesS 

Ai - A a = $(0, 0, a 3 ) - 0"d *(0, 0, a 3 ) > 0 (144) 
that is 

0(0, 0, agj (l4g) 
*(0, 0, B,) 

Then the equal sign, replacing the unequal sign in (144) 
and (145), gives the Wagner stress a M - The reason thfit , 
in spite of the doubtfulness of the method of making ar- 
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hitrary assumptions for several variables, the errors are, 
on the whole, .not. excessive ,- iA-pxobjably .due. _to .the faot 
that In (145)- accidentally, at least - the numerator G 
was made, a minimum through (143), oecauso It consists of 
only two quotas: 



| / EC* cp"*" dz and M G- J T <p |B ds 

o o- 

and from 0* it oan readily "bo shown that undor all pos- 
sible values it has tho smallest. 11 ' For angle and T 
sections, for example, 0=0. 

Severting to figure 7 and equation (131 ), applicable 
in the case of symmetry with the x axis: While the min- 
imum of ar^ is equal to Op, (equation 58), the require- 
ment of a minimum value for the numerator gives, with al- 
lowance for equation (15): 

(x A )* = ^ = x M (146) 

i.e., the Wagner stress: 

0 J T + H| E (o - G J T + XtfiE C* 

(<Ja ) _ ff . \ (147) 

J p + P X M J P 

This gives the relation: 



* p 8 

0 = 0 - -r*— (148) 
J x 



valid for symmetry with the x axis^ which can he employed 
in the numerical determination of 0 . 

11 formulate: 6 [ / w* a <U-1( / w* dP) a ] = 0. Then .the 

jf I 1 

introduction of equation (14), representing a purely geo- 
metrical relation, gives the equations (10) defining the 
center of shear. 
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On the other hand, a much more comprehensive view 
of the relation existing "between the Wagner stress <jjj 

and the actual twisting- failure atresa than is afforded 
from equations (136) to (138) is obtained from equation 
(134b) if fly is expressed with qp : 



(149) 



Observing that 



Xu — III — Xi 



(150) 



the conventional sine formula gives the following relation 



G J in + — ~ K G 
* lf_ 

J T + 7 x M x D 



(151) 



For the numerical calculation of a D , this equation 

is not oxpedient, as it would first require the calcula- 
tion of tho distance Xtj botwoen cpntroid and center of 

rotation from the quadratic equation (130 ). The signifi- 
cance of this relation lies in tho convenience of compari- 
son with equation (147) for dir. Eoplacing X]) with the 

R 

shear center distance xjj = — 2 in tho left-hand side of 



equation (130 ), it "becomes negative because , obsorving 
equation (148), we obtain tho following inequation: 



0* + 



TT 



a 



< 0 



(152) 



x valuos which in amount are greater than the roots of 
equation (130) make, on the other hand, the left-hand 
side positive, giving a picture as shown in figure 10. 
x H accordingly ■ always lies between the two roots x D , of 
which one is positive, the other negative (their product 
is given through - J p /F!) according to equation (130 ). 

To obtain the lower critical stress qj, it is necessary 
to insert the value having the same prefix as Xjj in equa 
tion (151 ). But this xj) value being always greater in 
amount than xy according to the foregoing, it follows 
that Op < ay. 
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For great Jij I the distances x D and xy and 
. consequently, . a D .and ct m .. differ very considerably, ac- 
cording to figure 10. The difference "becomes less as the 
wall thicknesses or fl lengths become less. But in the ex- 
treme case of J™ I =0, they disaopear for angle and T 

& "*" 

sections only, "because G =0 for these sections only, 
so that then the equal sign takes the place of the < 
sign for J T l B = o in (152). For all other (not double 

symmetrical or point symmetrical) sections, thero always 
romains a finite distance "between contor of shear and cen- 
ter of rotation- and consequently, a difference "between 
Ofc and Ojj (soo also equations (139) and (140 )) oven if 
8=0 or 1=0. 

7. SUPPLEMENTARY NOTE FOB. THE CASE IB" WHICH 
I and b ARE VERY SHALL 



For Tory small column length and wall thicknoas tho 
theory advanced hore fails "because then tho twist is ac- 
companied "by "buckling of tho walls and the aross-soctional 
shape is, as a result, altorcd. The twisting-failure 
atrcssos as computed with equations (58) and (99), respec- 
tively, would in that case "bo no more than a very rough 
approximation undor cortr.in circumstances, and a. mathemat- 
ical troatment of tho oxisting conditions is quite diffi- 
cult. For curved- section contor lino, for example, it 
would present a complicated shell problem Involving other 
than the usually allowed-for "tensile stresses" a, 

"flexural stresses" or linearly distributed over the wall 
thickness. _ Besides, the immediate effect of the shearing 
stresses t on the strains would in most cases be no lon- 
ger negligible, as they are in noway always small on ac- 
count of the short column length (shell length). They are, 
therefore, in general, no longer computable with (30) from 
the equilibrium conditions (21) and (120). 

For sections built up from a few straight pieces, the 
solution may perhaps be somewhat easier, but here also, 
the argument made for f , holds true. On the other hand, 
allowance for the bending stresses of linearly distrib- 
uted over the wall thicknoas is made easier, since it is 
possible to apply simple formulas for the change of oross- 
soctlonal shape at buckling. The simplest caso of this 
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kind is that of a cruciform section. Here, evon for very 
littlo wall thickness (in first approximation) the cross-* 
sectional shapo is proservod and the twisting of the cross 
sections about the controid (double symmetry) is identi- 
cal with "buckling of the walls after one single half-wave, 
the walls being visualized as platos pin-ended at ono 
side and free at the other. The "tensile stresses" 9, 
constant over tho wall thickness, disappear in this case 
because C = C = C* = 0. If b denotes tho leg length 
(width of plates), formula (48), as well as Wagner's con- 
formable formula (77) (S = Ml), gives: 

cr D = 1_£T = E if (153) 
33 Jp 2(1 + n) b a 

Tho twisting-f allure stress or buckling stress ctjj 

is accordingly unrelated to the column length. But this 
result io valid only if the length does not fall below a 
certain amount. Because, as tho column length becomes 
less, tho hitherto neglected "flexural strossos" 5 lin- 
early distributed over the wall thickness , become of 
ever— increasing significance and result hore in a rise of 
critical compressive stress. The completo Wagner formula 
(soo footnote, p. 32) allows for those bending stresses 

a by adding tho torm c odn » related to tho transvorse 
warping, to tho term C D £ U = C*, connectod with the cir- 
cumferential warping and tensile stresses flf , respective- 
ly. Eoro C tiu = 0 and C^ dn = (J) b 3 s 3 ; and (153) is 
replaced by the more exact values: 

-w- - C-crV^y + S ?) ' (154) 

This ocuation is, up to a factor 1 : (l - \i B ) in 
tho correction term. In agroomont with tho. approximate 
term originally obtained by Timoschenko (roforonce 5 ) on 
the basis of the Eitz nothod for buckling of platos under 
tho previously cited boundary conditions. It is seen that 

tho improvement through ^-^^ n for ^ > 10 is already in- 
significant (< 2 percent). 

According to Wagner, the same formula (154) would be 
used on an anglo section of leg length b, but this prop- 
erly would be permissible only when the axis of rotation 
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D paa8od through the point of intersection of the tiro 
logs* And this really nevor happens, ae shown at the . end. 
of-the preceding chapter (unless it is Intentionally ef- 
fected through guides, etc.)., while on the other hand, the 
error is so much less as wall thickness and length are 
loss. 12) ji» or greater ratio l/h, the porreotlon given 

through (jLjj) (^) *s meaningless, 'because then, proper- 
ly, it would require a correction toward the other side, 
which allows for tho fact that center of shoar and center 
of rotation arc not coincident. Even so, the orror is 
less than 11 porcont despite the great value Jx^j/ip* = 

0.61, as can ho readily proved, "because at fairly small 
ratio l/h, Euler buckling already takes plaoo in direc- 
tion of tho axis of symmetry ^at (iQ* Q^* * 1.07^. 

Similar conditions prevail on the^uneven leg-angle 
section and on the T section, where C likewise = 0. 
But for all other sections, allowance would have to he 
made for tho chango in cross-sectional shape induced. by 
buckling. Evory section therefore presents a new problem 
which is probably solvable for the simpler cases only. The 
actual critical stresses r.ro lower, rather, than higher, 
comparod to tho twisting-failure stresses computed accord- 
ing to (58) or (99). 

8. PRACTICAL DETERMINATION OF THE VALUES Ry , R x , and C 

The starting point for computing the warping moments 
Ry, R z and of tho warping strength C dependont on the 

cro ss-so ctional Bhapo only, is formed, as seen from equa- 
tions (37) to (39) by the unit warping w referred to tho 
centroid which, according to figure 1, can bo determined 
as aroas. Por sections built up from straight piecoa, the 
calculation becomes fairly simple booause the w aro in- 
flections, linear functions of the circumferential coordi- 
nate u, as well as tho coordinates x, y of any point 
of tho section center lino. Thon the integrals for Ry, 
R z , and C can ho computod exactly from. -the following for- 
mulas, in which i|f and (u denote any linear function of 

TO \ 

'In contrast wJLth the doubly symmetrical cruciform sec- 
tion, 5 and C disappear on the angle section only in 
the extreme case. 
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ti at interval Uk-i % u % u: 



^k 



\|/ (u)u) (u) du = — ^ [♦fc.^ (a^k-i + <" k ) + 

U k-1 + *k (2w k + CUk-i)] 



/ 



+ u> k (2^ k + i|r k-1 )] (155) 



\|> S (u) du = [^ k-1 3 + * k-1 * k + ¥ k B ] (156) 



Uk-i 

Those formulas, whose calculation is "best effected hy 
tabulation can, of course, he enployed also for ascertain- 
ing tho lnertial and centrifugal moments. They are also 
applicahle for partially or completely curved center lino 
if it is replaced hy a set of straight s.- 1 - 3 ' 

Quite often it will he of advantage to, apply, instead 
of the w referred to S other unit warping w referred 
to any favorahly located reference point 0, With the 
equation (fig. 2): 

w = w + x B y-y 8 3C+K (157) 
a simple calculus gives the following relations: 

By = x w dF = Ey + y s Jy - x s ( J xy + j v dl 1 ) (158) 

E x = J 7 » « = H 5 - s s J x + y s (j xy - x? dp) (159) 

13) ■ 

For the usual section forms, a tahulation of the numer- 
ical values for , Hy , C would he suitahle similarly' 

as for tho cro ss-Beotional aroa and the inertia moments. 



. 851 ' 43 



2y fl Ry (160) 



dJ (161), (162) 



The application of these formulas to a channel sec- 
tion will sorve as illustration (fig. 11), Owing to the 
symmetry with the x axis, it is Ry = 0. lor computing 
R s and C, it is "best to use the coordinate system z, y 

parallel to x,y that is placed through the intersection 
0 of section center lino and axis of symmetry. Observing 
tho prefixes for r^ and r^ as given in section 2, it 

is seen that, when passing from point 2' to point 2, r^ 
has the positive value a/ 2 for "both flanges and value 0 
for the web. The coordination of the unit warping w Q a 
0 to point 0, is followed according to equation (101) 

u 




o 



"by the unit warping antisymmetrical to the.- x axis! 
w a i = - ^jr, w^i =0, w 0 = 0, Wjl = 0, w a = + t£ (163) 
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C = J w» dP - j (J w dl) = 0 + x fl B 



f 

- 2x B y B J iy + y a a J y - 2x B % + 



Hereby 



Ry = / x w dF, R£ = / y w 
F F 

and, corresponding to (107), 

C m J w* dF -± ( /'w dF/ 
' F ' J 



whence 
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14) 

w dP = 0 



S 

Then equations (155) and (156) givo ( In simple fashion: 

B= = - 7 s a a h a (164) 
4 



and 



It further is: 



and 



C = + - s a 8 h 3 (165) 



x B = d : (a + 2h) (166) 



J * = 12 8 (a + 6t) (16?> 



The final rosults, according to equations (159) and 
(160) aro the oxprossions: 



ffx _ tjstjt . (168) 



p n a. - 8 t o- * s a a h 3 2a a + 15aD + 26h a , * 
C = 0 + x 8 J x - 2x a R-£ = -jg ("a-TlT)* < 169 > 

previously omployod in the oxanple of section 3. 



Translation "by J. Vanior, 
National Advisory Conmitteo 
for Aeronautics. 



14) 

Choosing the opposite sonso of rotation changes the 
profiac of r t , hut not of w. 
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Figure 3.- Determination 

of the twisting 
failure stresses tfp and 

2 3 

from o" s ,<f x ,p x by means of 
liohr's circle. 



e m 



Figure 6.- Channel sections With 

guide s ( enforced axis 
of rotation A) . 



Figs. 1,2,3,4,6 

Corresponding 
to: 

V-i-fP I P-Ff 

\du cosa-'iy 
(x,jfj \au sina-dx 



Figure 2.- Derivatives of 

■ the unit warping ■"*• E > 5 - 
x-y-System through T, - B , + r„ du 

centroid S; j 0 J 

x*-y*-Sy stem through 1 „.=„,.. + f r( . du 

center of shear M; J o 
x-y-System through any -\ _ _ P " 

reference point 0; j 0 J • 

x-y-System through actual! 

center of rotation D;J 




= S 0 + J ?i du 
o 



~( = rc + X s COS a + 5j sin a 
J r ( = COS a + Vn Sill a 

r i = rt + S s cos a + i/.s- sin a 
r t = rt" + * « cos « + sin a. 




Figure 4.- Twisting failure stress 
Oj) against length 1 and 

wall thickness s of channel sections 

with a/b = 1. 
Euler stresses (f x and ~dL according 
to Eq. (73), (74), twisting failure 
stress tfjj according to Eq.(58) by 
means of Eq. (73) ,(75) ,(76) , Wagner 
stress ct^ according to Eq.(78), 
Bleich stress ^igij,}. according to 
Eq.(79) . " ^ 
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Pigs. 5,7,8,10,11 
Figure 11.- Solution 

of and 
C on XT-channel section. 
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Figure 5.- Chart for determining the 
twisting failure stress 
from the cubic equation 

<u» — 3co>+pco — g=0. 
n a -j- ax + gy 



j» # # - V p es o as to js & 

Figure 7.- Twisting fail 

ure stress 
with enforced axis of 
rotation A for a channel 
section with a/b= 1, 
l/b= 24, s/b= 0.1(Fig.4) 
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Figure 8.- Curves for computing 

twisting failure stress 
o^j from Wagner's stress ojj and 
Euler's stress <f x for open sections 
with symmetry to the x axis. 
X|j/i*p for channel sections from 
Fig. 9; for other sections compute 
with Eq. ( 15) . curve for chan- 
nel sections in the extreme case 
(5j» <fy. For channel sections with 
a/b>D.73, where C <cr x ,it gives 
the maximum values y of o^/op. (Digits 
denote length-width ratio a/b). 
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Tigure 9.- x 2 */^, 2 and J^j /l^ 

values of channel flections. 
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